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Abstract 

We calculate the zero point energy of a massive scalar field in the back- 
ground of an infinitely thin spherical shell given by a potential of the 
delta-function type. We use the zeta functional regularization and express 
the regularized ground state energy in terms of the Jost function of the 
related scattering problem. Then we find the corresponding heat kernel 
coefficients and perform the renormalization, imposing the normalization 
condition that the ground state energy vanishes when the mass of the 
quantum field becomes large. Finally the ground state energy is calcu- 
lated numerically. Corresponding plots are given for different values of 
the strength of the background potential, for both attractive and repulsive 
potentials. 

PACS class 04.62 +v 



1 Introduction 

In the recent years much interest has turned to the study of the Casimir energy 
for spherical configurations. The issue concerns in particular the attempt to 
explain by means of QFT the puzzling phenomenon of sonoluminescence [Q], that 
is the emission of short intense pulses of light by collapsing bubbles of air in 
water. Up to now the hypothesis that the Casimir energy could play a role in the 
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photon emission has not been supported enough by a satisfactory QFT-model 
for the dielectric ball. To calculate the Casimir energy for such a ball could 
be of great importance to compare, in the statical approximation, the order of 
magnitude of this energy with the effectively observed light pulses . Many authors 
have already studied a perfectly conducting spherical shell in the vacuum of the 
electromagnetic field, we cite here only paper ||, other authors have investigated 
the conducting shell for massive fields [[|. In this paper we propose the more 
realistic model of a semi-trasparent shell. In calculating the ground state energy 
we use zeta regularisation techniques and Jost function approach. The heat 
kernel coefficients for a penetrable sphere have alreay been calculated in a very 
recent work Q . We carry out the complete calculations in order to find the exact 
numerical dependence of the ground state energy on the radius of the sphere. 



2 The model 

We want to study the ground state energy (GSE) of the scalar field ip(t, x) (to be 
quantized) in the background of a potential V(r) . We will consider the following 
field equation 

(□ + m 2 + V(r))<p(x) = 0, (1) 

where m is the mass of the field. The spherical shell is a geometrical object with 
radius R and a surface S, to whom it can be associated a classical energy in terms 
of classical parameters. The total energy of the system reads 




where V is the volume of the sphere, p is the pressure, a is the surface tension and 
F, k and h are other parameters. The classical energy is expressed in a general 
dimensionally suitable form which depends on powers of R, this definition is use- 
ful (as we will discuss later) for its renormalization. The quantum contribution 
in (2) is the traditional expression for the vacuum energy of a scalar field whose 
energy eigenvalues are u n . To render the eigenvalues of the energy discrete we 
take a finite quantization volume with radius L ^> 1. 

The classical shell is static and spherically symmetric. It is described by a po- 
tential 

V{r) = %S(r-R), (3) 

where a is the magnitude of the potential. The semi transparency of the boundary 
will be imposed later. The quantum contribution to the total energy is divergent, 
for the regularization we adopt a zeta function technique. We define a regularized 
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ground state energy 

^ = ^B A (n)+- 2 ) 1/2 -v s . (4) 

(n) 

where fx is an arbitrary mass parameter, s is the regularization parameter which 
we will put to zero after renormalization and A( n ) are the eigenvalues of the wave 
equation 

[-A + V{x)\<j> {n) {x) = \\ n) ^n){x). (5) 

Now we introduce a zeta function. The zeta function of the wave operator with 
potential V(r) as defined in (5) is 

CvOO=$>M+™ 2 )~ S > ( 6 ) 

in) 

then we can express the ground state energy in terms of the zeta function 

E, = \Cv{s-\)n 2s . (7) 

Since is valid the equalty 

1 r°° 

X " = W)L itt " le ~"' (8) 

we can write eq.(6) in the following form 

1 r°° 

fr(*)=£f£)/ dtt s -'e- x ^-™ 2 \ (9) 



(n) 



that is 



1 poo 

Cv(s) = — dt t-'e- mH e ~ X(n)t > ( 10 ) 

1 [S) J ° (n) 

(v(s) = f^l dtf- l e- m2t K{t). (11) 

The function K(t) represents the heat kernel. Now taking the asymptotic expan- 
sion of the heat kernel for t — > 

/ 1 \ 3/2 00 1 
K(t) = ^exp(-A? n) t)^° — J = 0,2.1- (12) 

(n) ^ ' i=0 

and making the substitution s — > s — 1/2 in eq.(ll), we get an expansion of 
in which it's easy to isolate all the divergent (pole) terms. This makes possible 
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to define a total divergent contribution to the ground state energy 

E - = ~64^U + ^~2) Ao ~24^ Al/2 
m 2 ( 1 , 4u 2 \ , m 

1 ^i + ln^-2W (13) 
s m 2 J 



32tt 2 



The Aj are the heat kernel coefficients. In the definition (13) we have included 
the fractionary terms A\ji and A3/2, which do not contain poles, to satisfy a 
renormalization condition (see below). Actually only the first three terms seem 
to contribute to the ultraviolet divergencies. In a smooth background potential 
their corresponding cofficients are well known: 



A = j d 3 x ; 

A 1 = - J d 3 xV(x) ; 

A 2 = X - J d 3 xV 2 (x). (14) 

The renormalized zero point energy can be defined as 

E r ; n = E V - E*». (15) 

To keep the total energy of the system unchanged we add the subtracted object 
to the classical energy. Then we have also a definition of a new classical energy 

^class E c i ass 

+ K m - ( 16 ) 

The transition from E c i ass to e c i ass consists in the renormalization of the classical 
parameters contained in eq.(2). Since the heat kernel coefficients are geometrical 
coefficients depending on the background (and in our case containing powers of 
R), to each classical parameter in the classical energy will correspond dimension- 
ally a term in E^ v , then we renormalize the parameters as follows 

m 4 (\ , 4/i 2 1 
p — > p - - + m 



m 

a — 



64vr 4 \s m 2 2 

3 



24tt 3 / 2 ' 

m 2 ( 1 , 4u 2 
F -> F + 7 - + In ' 



32tt 4 \s m 2 
m 



k ~" * + 16ir3/2 

h - h -M\ +i ^-^- (i7) 
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However in our particular case the Aq and A\/2 coefficients, corresponding respec- 
tively to p and a , will turn out to be zero, then only the last three parameters 
in (17) will undergo renormalization. These new parameters are to be inserted 
in the classical energy in (2). The old classical parameters are infinite quantities 
(i.e. unphysical) because they don't take into account the vacuum fluctuations. 
Before applying this renormalization procedure we must note that the ground 
state energy proposed in (15) has not yet a unique meaning. For the uniqueness 
of E^ en we impose the renormalization condition 

lim E r ; n = 0, (18) 

m— >oo r 

which physically means that for a field of infinite mass we have no quantum 
fluctuations. We fulfil such requirement by subtracting all the contributions in 
£div proportional to positive or zero powers of the mass. That is we subtract also 
"fractionary, terms up to the term resulting from the heat kernel coefficient A-i- 
The remaining part, containing negative powers of m, will go to zero for m — > oo. 
This is our renormalizaton scheme, now we turn to the calculations. 



3 Representation of the ground state energy in 
terms of the Jost function 

We use an approach appeared for the first time in 0. In the background of a 
spherical potential we have a radial Schrodinger equation 

d 2 1(1 + 1) 



dr 2 r 2 



V(r) + \ n i Mr) = 0, (19) 



where I is the angular momentum. In the general scattering theory with a conti- 
nous spectrum p we have the "regular solution,, [0 defined as 

<t> P ,i{r) r ~jl(pr), (20) 

where ji(pr) is the Riccati Bessel function. The asymptotics of the regular solu- 
tion is expressed in terms of the Jost function fi(p) 

Mr) ™ \ (fi( P )K(pr) - ft(p)h+(pr)) , (21) 

where hf(pr) are the Riccati-Hankel functions. Now we examine the field at 
the boundary of our quantization volume. The potential has a compact support 
then, at the boundary, expression (21) becomes an exact equation, which can be 
considered as an equation for the eigenvalues p = A n ^. Now taking for instance 
Dirichlet boundary conditions: <f> p j(L) = 0, we will have 

fi(p)k(pL)-f! l {p)ht(pL))=0. (22) 
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We see that eq.(22) is satisfied for p = Xi jTl , then we can write the sum in (4) as 
a contour integral using Cauchy theorem 

^ln (fiipfcipL) - ff(p)k+(pL)) , (23) 

where the contour 7 encloses all the solutions of eq.(22) on the positive real 
p axe and also the bound state solutions in the limit L — > 00, which lie on 
the imaginary axes . We further simplfy eq.(23) by separating the contour in 
two pieces 71 and 72 and developing the Hankel functions for large L. Then it 
is possible to recognize in the integrand a term ipL which corresponds to the 
Minkowski Raum contribution. This term can be dropped. Now we shift the two 
contours 71 and 72 to the imaginary axe and substitute p — > ik, so that we get 



_ < ~^ 1 /*oo O 

^ = -^/x 2s ^(/ + l/2) / dk[k 2 -m 2 ] 1/2 - s ^ In fi(ik) 
71 1=0 ^ m 



(24) 



Since at the end our quantization volume will go to infinity ( L — > 00 ) this equa- 
tion will be independent from the boundary condition chosen for the quantization 
volume. Eq.(24) is a very general and usefull representation of the ground state 
energy, where all the information about the background potential is contained in 
fi(ik). 

In order to perform the substitution s = and the subtraction proposed in (15) 
we split eq.(24) in to two suitable parts, one of which is divercence free. We 
obtain this by adding and subtracting the leading uniform asymptotics of the 
integrand in eq. (24) (for more details on this procedure see for instance ||). We 
define 

E r ™ = E f + E as , (25) 
Ef = -^V E( Z + \) £ dk l k " - - Infrm] (26) 

and 

Eos = -^/i 2s £(J + \) J°° dk[k 2 - m 2 } 1 / 2 -^ In - Ef, (27) 

where f[ S {ik) is the asymptotics of the modified Jost function which we will take 
up the third order in v = I + 1/2. As the reader can check, the quantity E { 
introduced in (15) remains unalterated. 
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Now we need the modified Jost function corresponding to our scattering prob- 
lem, so we turn to study the potential. 
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4 Jost function of the £-shell 

The initial field equation 

(□ + m 2 + V{r))ip{r) = (28) 

valid for — oo < r < oo can be divided into two parts: an equation for the free 
field 

at r^R — ► (□ + m 2 )y?(r) = (29) 

and an equation for the field on the shell, which includes the required trans- 
parency conditions 

a t r = iJf fl + >-^ B - >= ^ (30) 
[ — > continuous 

For the delta potential the regular solution is 

Mr) = ji(kR)Q(R - r) + ~ ~ flight (kR)) 0(r - R) (31) 

for the field inside and outside the radius R. As above ji(kR) is the Riccati-Bessel 
function and hf(kR) are the Riccati-Hankel functions. Combining eq.(30) with 
eq.(31) we get 

Mr) = Z -(f l (k)hr(kR)-ft(k)hr(kR)), 
^jiAr) = k l -[Mk)kf(kR)-ft(k)hf+(kR)). (32) 



We solve for fi(k), keeping in mind that the Wronskian determinant of hf is 2i, 
then we find 

/,(*) = i (-2z(-l) + 2l ^UkR)ht(kR)) (33) 

or 

o 

f l (k) = l + — ]l (kR)hl(kR) (34) 

and for the modified Jost function we get 

f v (ik) = 1 + aI„(kR)K u (kR), (35) 

which is in terms of the modified Bessel functions I u and K v where v — I + 1/2 . 
We need also the asymptotics of the Jost function: f® s (ik), or more exactly the 
logarithmus of f^ s (ik). The expansion of the product of the two Bessel functions 
in (35), for k and v equally large, is easily obtained with the help of |§. Then 
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we find the needed asymptotics as a sum of negative powers of v with coefficients 
Xj >n depending on a. The asymptotics up to the third order reads 

at a 2 t 2 at 3 a 3 t 3 3a t 5 5a t 7 /l\ 4 
2 V ~ Y^ 2 + 16^3 + ^7^3 ~ + ^6 ^ + UJ ( ^ 



Here as above is v — I + 1/2 and t = l/^/l + Now the ground state energy 

j^ren _ _|_ can explicit calculated. 

Let us first begin with an analytical simplification of E as . We transform the sum 
over / in an integral with the help of the formula 



=o 



valid for a each continuous function F{y). In our case is 

F{v) = J dkv\k 2 - m 2 } 1 ' 2 - 3 ^ In f?(ik). 
Therefore E as is split into two addends 



7T 

and 



oo 



E^ = -—y?-l duF(u) (38) 



o 



EjS = _c^sns f°° du (F(iu)-F(-iu)) 

ix J n 1 + e 2 ™ i v ' 



For the integrations we use the formula 



Jo L 1 J dku» 2 r(|) 

(40) 

First we calculate E^: 

E<$ = -?^V/ Avv \ dk[k 2 -mY' 2 -°?-\nf?(ik)) 
it J J m dk 

m l - 2s ii 2s \ ^ (mi?) 2 "™ r(| - s)r(l + *=?)T{s + nz2) 

2^ A .'-"" 



7T 2 r(|) 



Here, inserting the coefficients of (36) and expanding up to the first order in s all 
the terms which depend on the renormalization parameter we get 

Ks ~ 96nR L m 2 2 



Rm 2 a ( 1 , 4u 2 
' +ln ' 



87T V S m2 

ma 2 
+ ^6~' 

the terms proportional to m 2 and m° contribute to the divergencies of the energy. 
They are used to calculate the heat-kernel coefficients in (13) and will disappear 
after the subtraction of E^ v . The term proportional to m corresponds to the 
A3/2 term of the heat kernel expansion; although this term seems to generate no 
divergency it will be as well subtracted as mentioned in the second section of this 

paper, because of our initial renormalization condition . 

(2) 

Now we calculate E as '■ 

' m i-2 SfX 2s\ / 3 \ r ( s + / rR dv v ^-«2 cos[?r(j -n)\ 



(=^) r 



/ ^ ": 2 .r M S: s)1 1 . (4i) 



then 



32vr 2 Vs m 2 J 2nR U ; 
h<2,{Rm) 



a 2 



8R 2 m 
2a 3 + 3a , /n . 

h±\Rm) 



8nR 

where h n (x) are the following functions coming from the ^-integrations in eq.(41) 

f°° v v 2 
hi(x) = / dv — In 1 1 -I ; 

V ; J 1 + e 2 ™ 1 x 2 ' ' 

f x v 1 

h 2 (x) — I dv 



1 + e 2nv 1 



f°° . z/ 2 . / v 2 \ ' 

poo 

h±(x) = / dv\n\l — 



v 2 , fl ( v 2 ' ' 



x 2 ' \v V 1 + e 2 ™ 
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*.(,) = I d „ H l-_H^_^ j . (43) 

(2) 1 

We see that also E as have a pole of the form = for j = n = 1. This pole will 
contribute to the heat kernel coefficient A 2 . 

Now, since Ef contains no poles, we are able to write down the complete heat 
kernel coefficients Aj, up to the order j < 2: 



A = 0; A 1/2 = 0; 

2 ct 3 

A 1 = -AirRa; A 3/2 = vr 3 / 2 a 2 ; A 2 = -g^r— ■ 

This coefficients are the same as in paper [f|]. After performing the subtraction 

(Eg + Eg) - E^ 

Ei 1 } cancels completely and only E^} (whitout its divergent portion) will con- 
tribute to the total energy, then we have 

oi 

E as | s=o = -tt - ^ hi(Rm) 

Z7lK 



a 2 



h 2 (Rm) 



8R 2 m 

+ 2a3 + 3a h 3 (Rm) 
ISttR k ' 
a 

-— h 4 (Rm) 

+T— 75 h(Rm) . (44) 



487T.R 

5 Asymptotics of E as 

As a test for the plots that we're going to do, we check analitically the behaviour 
of E as for small and for large values of R. To do so we must find the corresponding 
asymptotics of the h n (x) functions. 
For R — » we find: 

lim hi(Rm) ~ Infirm) ■ h Ci ; 

hm/i 2 (Rm) ~ R 2 m 2 + C 2 ; 

lim h 3 (Rm) ~ -21n(Rm) ■ ( — j + C 3 ; 
lim h A (Rm) ~ -2 ln(i?m) • (-1) + C 4 ; 

R — >0 

lim /i 5 (ifcn) ~ -2 ln(ifcn) ■ (-4) + C 5 , (45) 

— >0 
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where the C n 's are numbers resulting from the //-integrations. Then we have 

]frE as ~--^--??-(h 1 ±-C 3 )+o(±-)\ (46) 
fl-o 16nR 24nR \ Rm 7 \R J 

where C = (-8C1 + C 3 - C A + C 5 /3) ~ 0.224 and C 3 ~ 1.96. 
For R — > 00 we have 

1 7 



then we find 



lim hARm) ~ 

u ; m 2 /? 2 1920 

lim /i 2 (i?m) ~ — ; 

lim h 3 (Rm) ~ — 1 ; 

7 

& W ~ (47) 



la 2 /O 2 



lim £ as + O — . (48) 

384 mR? \R' 



6 Numerical results 

To study numerically the behaviour of E as we write it in the form 

a 

E as = qARm) 

a 2 

+^(^») 

+-^^ 3 (i2m), (49) 
where the three functions g n (x) are given by 

#1(2;) = -8/11(2;) + /i 3 (x) - 2/14(2;) + -/i 5 (z) ; 

g 3 (x) = \h 3 {x). (50) 

As it's clear from (49), for small values of a E as will behave like function gi(x). 

For large values of a, E as will behave like function ^3(2;). 

The plots of the three g n (x) functions are shown below in Fig. 1 
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For the complete quantum energy we still need the Ef contribution. In the 
expression (26) for Ef, after putting s = 0, we integrate by parts and we obtain 

1 00 1 f°° k 

E f = -^ l+ oH m i Q*Mik)-hif?(ik))dk (5i) 

this quantity can't be further analytically simplified. Below we show (Fig. 2) a 
plot of R • Ef as function of it! for a — 1. 
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For the total ground state energy as a function of the radius of the shell we 
get the curves shown below (Fig. 3-5) for different values of the magnitude of the 
potential a. 




-0.6' 



Figure 3: The renormalized vacuum energy E^ en (R) for positive and negative 
values of the potential. 



0.002 - 
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Figure 4: The renormalized vacuum energy E^ en (R) for a=0.3 

7 Conclusion 

We have obtained an explicit expression for the renormalized ground state energy 
of a scalar massive field in the background of a semi-transparent shell. This 
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Figure 5: The renormalized vacuum energy E^ en (R) for a=10 



expression is given by the sum of (49) and (51) and it depends only on the two 
parameters of the classical system, namely the radius and the magnitude of the 
potential of the spherical shell. The plots of E™ n as a function of the radius show 
that for repulsive potentials the renormalized ground state energy is positive only 
in some limited intervals of the radius and only when the potential is smaller than 
1. For a potential larger than 1 the energy is always negative. This is the most 
striking conclusion of our work. For very large values of a the shell should be no 
more transparent and the problem should formally become a Dirichlet boundary 
condition problem. One could check this in the equation (35) for the Jost function: 
here inserting a large a the addend 1 becomes negligible, then we would just have 
the product of the two modified Bessel I and K functions; such a Jost function is 
exactly the one obtained in || for a perfectly reflecting spherical shell ( Dirichlet 
boundary conditions). In that paper the ground state energy is simply the sum 
of the energy inside and outside the shell. Then formally we should have: 



Now it is shown in || that the "mirror,, configuration has always positive 
ground state energy for repulsive potentials. This fights with our plots which 
show an opposite sign. Furthermore the Ai coefficient in paper || is zero. In our 
work A2 remains a non zero coefficient also in the limit a — > +00. 
Nevertheless for flat parallel boundaries the formal transition hypothized in the 
limit above is actually fulfilled as shown in paper 0. At least at the present 
stage of the investigation it is difficult to interpret these results. We must stress 
however that a perfectly reflecting boundary is not a realistic object and that a 
semi-tansparent sphere is a better physical model. 
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